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0 0۴ 02880 متقصتان 01 ۳۲1۲6۵۱۵ 10 .(2024) بک رتلم۱۱620۵0212 عک و.۲ متطفاه2ظ 
0۷۵۵۵۰ صمم مک ۵۶ طنامو :نگ عققن)) )صقاظ مفقا۵ظ ۵۶ وععه)9 ل2منع10مععظظ 8ه2 
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۱۱۱۱۹۹ 


عافصته ۶ معصه)0۲مص1۳ معط مهن .اعزم۵و صححصاط ۵۶ فصهعصمع صه)1۳01۲ دمح فص ۵۶ عصم ود واتتتاععو ۳۵0۵0 
21 ۵۶ )20۷76 م1 ,لقتاصمووم ود همه فتطا صا فافهععع10 مافصصئاه مصما-عمما مصانمنطامه رتماععی آهعتتلنه‌توه مطا 10۲ 
عافصتله 0۶ امجمصصا عط) رنه ما فعمطل‌نهعوع2 عم هاماتوومم 11 فص فقط (۴0) واممصد تمطمنوع هه (00) متقحصتاه 
مج 100106 مایق متتمطامومصصه صا 286عععص1 فطل «القلهعهروع رمعصمط متفصتاه ۲۵00)60ظ .امعم لهعتلنه‌تتهه وه 
موه وه وا0عه مرم۵۷و 2۷۵ از بالحگصنه: ۵1 مصمانهم مه تصنامصصح فطا مد ومعصقط موه ۷۵1 وه ,مها 
۵ م3 طا ءعصفل بح رعع۵10عظ1 .فممتوع 0زعه-تصصمه مه 2710 صا «الهتمعووه ,انعر 1۵00 مه ممتامنالم۳م 
,۲صمنوعوصهن رصح ۷۵۲۱۵ 6 ۵۶ مادم اصممگن صا الم تمتنملنم‌تنيه ۵۶ ممتامن0۲0 مطا فامعه تمعن رمتقصتاه 
ععصقطه متفه ۵1 فامع]له مطا ماهمتادع1۳۷ م1 6000۵60 ۵۵ له ونط 1 .۱9۵۲10 فاص بماتتناهعع 1000 ۵۶ یاوعد طا 
۰ 1۵۲۳۱۵2۴ 0۴ فطمنوع۲ لهمزو۵ماناه لمع آهعزم۲۵ طا صقان ملقامم ۵۶ ومعهاه لهع1ع10مصمدصام ۵۶ طاعجع1 0ج ۱۷1610 


1۷]2)6۳۵216 20 ۹ 


,(2005 60 1989) وعدوع۷ 17 0۶ 0۵04 )0209و 2 طاذ فصمتنهاه وع0امصا طم‌توموع۳ فنطا ظ همه تناو فط1 
0۳601012408 ر22) عاجهم‌جطهم مصیتمصتصتدطر لمح تنم رولبمه فنط م1 مصعزمصه]۱۷ فصه رزممصطفک باگ0ت1ل نت101 
06 10۲ 2۵0 ۷۷۵۲۵ م1هعو محصتا «انعل ه وه مط وه 4ص رلععصو مصه معهعع7ه ,اتلتصنط 10120۷ و2۷ 
۶6 136 مها ۱۵۵0 ۷28 [06مصظ صمتهانهت قمتمجعع ۲۲2قطصهن فطا ۵۶ آبامتان فطع حصم علهعوممنمه أهتاصمصمهرم 
۶ ]000۵ عط 61و مدمه ما 60و وه عم قملهعومتصه او مظ1 ,۲۵۵8.5 208 ,۲0۵4.5 ,۲۵۴۵2۰6 808821109 
عطا ۶ انلتله۷ عطا مصتصیعاع0 ما مه ممتاه‌طاتلم 10۲ مه 116010 منهجععج مه عع0ع۵ م1 .061مقظ صمتاهانمته تقعصهع فطا 
تمعن هویم لمتانهم‌تنوه ما ۵۶ فامعزمزن طمندعوع۲ مت ماه رفصمتاتجم لفتاصعامم ما عم ]۱۷۷۵۳۵۵ 
,0216 عصتاصهای ععلبامصا ممتاقصم‌کص فنط1.. ,عنام ملفهه 2 فه 4عفا 24 0ع)عع۱آهع 76۲۵ 2010-2011 مصتتتال 
اممزع0۱0وطه مه وعممموممصه مه رممته۳2 10 ما عمصمعممه رععظمعتمصه ها عصتاصهان صم ومعفاه 2ع1عم1مصمدام 
٩017 ۰‏ عظ) ما کصفام متقاهم عطا ۵۶ تام 


جمتکمی وز(1۱ 20 وااباعع1 
عاصم200606 0۲0۷۱060 ۷۱۵۱4 ماهامم هه وععهاد لهمنع10مصعطام تم ۷102610 اع0مصه ۷۷۵۳۵۵ ۵۶ و)لناوع؟ مط 1 


جمتااطاتاظ عصمصصصصمت 0۵2۷۵ تفص 0ماباطاتتاوتل مامتاته وقععععه صممرم صنق و فتط1 ,(ع)مطانه 186 ۵2023 
مجح جمتاباماتماعنل مصمتامامم0ه رعوصتحطه بععه فاصم طمتطه ب(ر40 ۳۲۷ 00) عوصمعنا اقممتامجم‌م 4.0 ۱( 0( 
(۲)6مطانج امصتعتره معط مه اتمه ماقتاممتمم2 2۷7۵ ۲۷۵ قق عع10 قه ماقم ۵ حصباتلعصه رصح صا ممتام۲۵0۲۵0 جج 
۰ 136 200 
6 تا و و( ِ 


۴ نشریه بوم شناسی کشاورزی. جلد ۱۶ شمارها بهار ۱۴۰۳ 


۵۵6 ]0212006 ل۱۵)۵0۲0102162 ماقلتادله ۵ اع0مصهر تاه عطا ۵ طتلتطاحه م1 .ععاصاهتیه۷ طاهها ۶۵۲ ومتفطتادم 
0ص صباصصت‌رمصه 0ماهلباصته م1 ,مماهاامتمعن عمتاهلیصته ما پرمه‌تتهعه فوع1 09۵0و اع0ممظ مه رته097۵۷ظ رلعصتومع 
عطا 0 عمجم 1062860 ومتتهجهعو نک تمل‌صا (2011-2100) ۵۵۲100 ۱62 فط) 0۲ وععتاهتعم‌جصع . حصیتصتصتاظ 
معط لاه طا منهج عبط صا معصقط للذ۳ موله فقط ممت)ه)ز۵۵01۵ظ .ومنتمصهعو ظ معط 2 ما منم متام تعوها۵0 
46ص امه متقصصتاه 8۷۲2صدن) فطا ۶ فقلناوع مط عصتصتطصم عم ممتصتامع0 ه سم لصه. ومتت2طعع8 
ممتاهان له همع ۷۷۵۲۵9 عظ6 ۳168 2071-2100 280 2041-2070 ,2011-2041 0611005 عاقصصتان مععط صاً 868821105 
کصمتانلجمه افتاصمامج ممصتا ۷1۵10 تمطابت مققامم 200 ومعهاه لهعلومامصمدام اصعععلنل ۵ طاعجع۱ مط هط 0760و 2۵061 
0 عبل روع‌هانت مطا م۵ جمتاهتتممظ آهعلوما0ملو ردان مط ما ممتاقصتصصهع 0۶ معماه 1۳ .قهعته «لتتاد مطا مد ۲۵0860 ۵۵ ۷۷7111 
ملقام 0۶ تلادع م1 .ععهاه فنط مد 0عتتتاهعن 0۵100 مطا 0۶ طاعصعا مطا صا ممتامنا0ع۲ اومهعزع معط رحمتاهسل بتمعجما واز 
امعطونط معط مفالنافع: عطا ما عصتل؟0ععض .ققعته لتقاه مط ص ۱۷1۵10 مد مومع 2 7۵0«مطو ممت)هانصته ۷۱۵10 بانط 
صمتاعصیظ تقانم‌صهاع ۵۶ ععماصهعتع م1 .(2071-2100) ۵۵100 70ظ) عط صر 0عتموهان قه7 عمصقصصم‌گتهم صا طمتام۲۵0 
فطل صا مفهعععصا مج 0760و فافع مطا رلحتعصعع ص .4.30 مضه 2 مه هه صم ۵0موهان ۸۵۵ طمتام۲۵0 
6 .6۲1008 معط مه مطا تمعن فتتوامی وععتعع0 0.9 00 0.37 ۵۶ فععههم‌جصع] میمصت تمصظ لصه مصتحصتصتمه 2۷6۲۵۵6 
0 06062860 فقط ب«اتسطممظ لهمتوماملورطان ما عمصمع‌عمصه مصم اصقام متقاوم معط ۵۶ وععهاه لهعلومامممصام مط ۵۶ طاعمع[ 
6 ۶ ومع( عمط ما متاع۲60 م1 .وما۲2وم‌مصعا تمهت نهص 4صه رتاحتصتظ ص2مص فطل مد عقهمتمصز فص ۱۳ ۵2۲21161 
م۵2 مطا 0صه رورم 1.7 ۲۵ 0.5 ممع۵ و1 ماد عصلیه110۳ مطا رورم 1.9 10 0.2 ههع0۳7ع0 15 عههاه عممتافصتمع 
0 ۲612160 و1 فعههاه لهعنع۱0مصعطام مطا مصتس ممتامنا0ع۲ ]ومامهتع م1 .واعل 3.5 6۵ 1.4 ههع0۳۷ع0 و1 تفع م6 مههاه 

۰ 3.5 0۶ ب«اتتتاقط م۲ 0226و عمصمعتمصه مطا 


صمتعصل(10) 


6 0۱ معصفطل مافصتاه 0۶ فاهعله ۵22076 مطا مد عفهعتمصا مج مهعت4ص رتتاه فتطا ۵۶ فالناوع: مط بلقته‌صهع 18 
,60162566 ۲۵11211 عظ) صرح معوزع معبتاجهم‌ججه مطا مه صهزمصها مه رزهمصطفک راگ0ع1! ۵۶ فعمتاد)ه معط ما کصفام متقاهم 
0 معمصدفل ۵۶ مامح مطا روقهععع0 فعماهاهم ۵۶ ۷۱۵۱0 مطا مضه ممفوعو عصا مه فطا ۵۶ طاعصعا مط رتالحساهه 280 
ععصنن رماههم‌جصما لمتاصصح مههه7ه صا وفععتعصا فاتوامن 062۲68 مصم ۶۵۲ ,مهن عفهعتمص1 71 وعتمصصتله ممتتان] 
۰ 0 16076266 ۷۷111 ۱71610 


۶ و9۵۵6 ۲۵۵۸۵۱021621 ,۱۵061 ۷۷۵۲۵۹۲ رمله۱۷۱۵۵5۵ أهتاممه۱۵۵ راع۱۷۵0 ممتامانهتتن لهعصون 1۵۲۱۷۵۲۵۵۰ 
۵ 1۵۱۱2 ۵۶ انا رم0)2ظ 


نشربه بوم شناسی کشاورزی 


7 ۵2۷۰۵۱۵۵ :۲۲0126۵0226 
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جلد ۱1 شماره ۱. بهار ۱۶۰۳ ص ۱۳-۸۰ 


اثرات تغییر اقلیم بر عملکرد و مراحل فنولوژیکی گیاه سیب زمینی 


(مطالعه موردی: حنوب استان کرمان) 


۲ باه ۱ نو ۳ اف 
تاریخ دریافت: ۱۴۰۱/۰۳/۱۵ 
تاریخ بازنگری: ۱۴۰۱/۱۰/۱۹ 


تاریخ پذیرش: ۱۴۰۱/۱۲/۲۷ 


چکیده 


سیب‌زمینی از محصولات مهم زراعی است که به‌عنوان یک محصول استراتژیک در تأمین امنیت غذایی در آینده شناخته شده است. با توجه به 
پدیده گرمایش جهانی و تحت تآثیر قرار دادن تولیدات محصولات کشاورزی» ارزیابی اثرات تغییر اقلیم بر تولید سیب‌زمینی در کشور ضروری می‌باشد. 
بدین منظور از داده‌های مدل اقلیمی 20۳81۷۲2 تحت سناریوهای انتشار 10۳2.6 ,10۳4.5 و 50۳8.5 در مناطق گرمسیری و 
نیمه گرمسیری استان کرمان» استفاده گردید. از مدل آماری ]81751۷7 برای ریزمقیاس نمایی خروجی مدل 20۳51۷2 و از مدل ۱۷۷۵۳۵۹۲ 
به‌منظور شبیه‌سازی عملکرد غده و مراحل فنولوژیکی سیب‌زمینی استفاده شد. در این مطالعه» مدل ۷۷۵۳۵۹۲ با استفاده از داده‌های مزرعه‌ای سال- 
های ٩۰۵۸۹‏ و ٩۲‏ کالیبراسیون و با داده‌های سال‌های ۱۳۹۱ و ۱۳۹۳ تعیین اعتبار گردید. نتایج به‌دست آمده از شاخص‌های آماری» نشان از دقت بالای 
مدل ٩151۷]‏ و مدل ۷۷۵۳0۵5 و مطابقت نتایج کالبیراسیون و تعیین اعتبار با داده‌های دیدبانی بود. میانگین مربعات خطا در دو مدل 129۷ و 
۲ در مراحل کالبیراسیون و صحت‌سنجی به‌ترتیب ۰/۸ و ۷/۴ درصد تجاوز نکرده و در محدوده عالی قرار داشست. ۸۷ درصد مراحل 
فنولوژیک و بالای ٩۰‏ درصد عملکرد غده سیب‌زمینی» در اطراف خط استاندار (97-) مراحل فنولوژیک و عملکرد دیده‌بانی شده قرار داشت. نتایج 
شبیه‌سازی در مناطق مورد مطالعه تحت سناریوهای )1 نشان‌دهنده افزايش دمای کمینه در محدوده ۰,۲ تا ۱,۲ درجه سلسیوس, دمای بیشینه در 
محدوده ۰,۲ تا ۱,۲ درجه سلسیوس و کاهش بارش در محدوده ۶/۴- تا ۳۷/۶- میلیمتر خواهد بود. نتایج حاصل از شبیه‌سازی عملکرد غده و مراحل 
فنولوژیکی سیب‌زمینی در دوره آینده توسط مدل ۰3۷70۳057۳ بیانگر کاهش عملکرد سیب‌زمینی به‌میزان ۰/۸- تا ۶/۳- درصد بود. میانگین حداکثر 
کاهش طول رشد با 7۳/۵ روز در سناریوی 56۳8.5 مربوط به مرحله رسیدگی فیزیولوژیکی (برداشت) بود. همچنین نتایج نشان داد که به‌ازای هر یک 
درجه سلسیوس افزایش میانگین دما سالانه. عملکرد غده به‌میزان ۴/۹۹ درصد کاهش خواهد داشت. 


واژه‌های کلیدی: ریزمقیاس نمایی» مدل گردش عمومی» مدل ۷۷0۳051 مراحل فنولوژیکی سیب‌زمینی 


مقدمه 
پدیده تغییر آب و هوا به‌عنوان یکی از چالش‌های مهم قرن حاضر 
در سطح دنیا شناخته شده است. تغییرات آب و هوایی 9 عواقب و آنار 


۱- دکتری» گروه جغرافیاء دانشکده ادییات» دانشگاه حکیم سبزواری» سبزوار ایران. 

۳ استادیا گروه جغرافیه» دانشکده ادبیات و علوم انسانی» دانشگاه جیرفت» جیرفت» 

(#:- نویسنده مسئول: ۴ ۵ 1۰۱۵2۳۵/212۵11 :اتفصط) 
6 ۱ و( 


آن از ابعاد اجتماعی-اقتصادی» زیست محیطیء سیاسی و امنیتی 
دارای اهمیت خاص می‌باشد و در ابعاد ملی» منطقه‌ای و بین‌المللی 
باید مورد توجه قرار گیرد. بخش کشاورزی به‌سبب تعاملات گسترده و 
مستقیم با محیطء بیشترین تأثیر را از پدیده تغییر اقلیم می‌پذیرد؛ و به 
تغییرات آب و هوایی و پدیده‌های خاص مانند خشکسالی» سیلاب‌ها و 
طوفان‌های شدید حساس است. ممکن است افزايش دما نقش مثبتی 
در تولید غذا داشته باشد» اما افزايش پتانسیل خشکسالی. سیل و امواج 
گرمایی. چالش کشاورزان را افزايش خواهد داد ( ,تحصعطعمطتصوظ 


۶ نشریه بوم شناسی کشاورزی. جلد ۶ شماره!. بهار ۱۴۰۳ 


2 در صورتی که اقدامات مناسبی در جهت کاهش اثرات تغییر 
اقلیم و سازگاری با آن صورت نگیرد این پدیده می‌تواند بخش عمده- 
ای از جهان» از جمله کشاورزی و منابع آب را تحت تأثیر قرار دهد 
(2014 وله ۶ تصعمتنصه؟ :ع12). بنابراین» هر گونه تغییر در 
وضعیت آب و هوای جهان به‌طور مستقیم در تولید محصولات 
کشاورزی مناطق مختلف دنیا و در نتیجه بر موضوع تأمين امنیت 

امنیت غذایی یکی از مهم‌ترین دغدغه‌های جامعه بشری است و 
پیش بینی می‌شود با رشد جمعیت برای تأمین مواد غذایی مورد نیاز 
جامعه جهانی در سال ۲۰۵۰ میلادی, تولید مواد غذایی حدود ۷۰ 
درصد افزایش یابد (2014 ,.له 6 01صمحع1(21). تغییرات پیش بینی 
شده اقلیمی به‌ویژه افزايش میزان دی اکسید کربن اتمسفر درجه 
حرارت و نیز تغییر در میزان و الگوهای بارندگی اثرات شدیدی بر 
خشکت و نیمه شنک کواهد داشت: آقلیم‌های خشک و تیمههعتک از 
یک طرف. بهدلیل ساختار اکولوژیکی خاص خود بیش از سایر قیم‌ا 
نسبت به تغییرات محیطی حساس بوده و آسیب‌پذیری بیشتری دارند. 
از طرف دیگر همانطور که آشکار است از مفاد بارز و آشکار پدیده 
تغییر اقلیم در بخش کشاورزی, تأثیر در اکوسیستم‌های طبیعی و 
فعالیت‌های کشاورزی می‌باشد (2019 ,22120). بنابراین» به نظر 
می‌رسد که وقوع تغیرات احتمالیاقلیمی در این مناطق اشرات قاببل 
ملاحظه‌ای بر سیستم‌های تولید کشاورزی به همراه داشته باشد 
(2016 مصقحصفاتنا )۱ 

علی رغم تلاش‌های بسیار زیاد اقلیم شناسان» هنوز عدم 
قطعیت‌های قابل توجهی در مورد اثرات بالقوه تغییر اقلیم در بخش 
کشاورزی وجود دارد (2012 ,.21 66 00206)). با این وجود» تغییر 
اقلیم بر تغییر فنولوژی و فیزیولوژی گیاه» بیلان آب در خاک تبخیر و 
تعرق و بارندگی موثر بر رشدء عملکرد و نیز یی گیاهان تأثیر درد 
(2016 ,طوت‌صه۲2 عک تامطه]گ2 12721). بر اساس تحقیقات 
صورت گرفته افزایش غلظت دی اکسید کربن اتمسفر سبب افزایش 
عملکرد گیاهان 3 مانند گندم. برنج و سویا شده و سبب افزایش 
(2007 ,له اه تعطم۳1). 

سیب‌زمینی (006200) گیاهی است یک ساله از تیره 5012026626 
با نام علمی (م 1۲۵5۲ 501070:۱) که اندام اقتصادی آن ساقه 


زیرزمینی می‌باشد (2008 ,6[600017). سیب‌زمینی 5۵1077) 
(م] :۵۲05 پس از گندم (م] 7111600 برنج 02۵) 
(0/۲۵.۲؟ و ذرت (. 16۲5 24), چهارمین محصول در سبد 
غذایی جوامع بشری است (2012 ۳۸۵0 که نشان از اهمیت نقش 
آن در تأمین پروتئین و نیازهای غذایی مردم دارد. سازمان خواروبار و 
کشاورزی-جهانی» سیب‌زمینی را یک محصول تامین کننده امنیست 
غذایی آینده در دنیا معرفی نموده و کشور ایران با داشتن اقلیمی 
متنوع و با تولید سالانه پنج میلیون تن سیب‌زمینی در رتبه سیزدهم 
تولید این محضول در جهان قرار دارد [2018 ,40 سیب‌زمیشی از 
نظر توازن پروتئین در غده‌هاء دارا بودن اسیدهای آمینه مهم سازنده 
پروتئین» ویتامین‌ها و مواد معدنی در تغذیه دارای آهمیت ویژه می- 
باشد. برخی از تغییرات رخ داده در دوره رشد گیاهی در مقابل تغییر 
شرایط اقلیمی» در مقیاس محلی و منطقه‌ای قابل مشاهده است. 

طبق مطالعات و تحقیقات انجام شده تأثیر افزايش غلظت دی 
اکسید کربن اتمسفر و تغییرات آب و هوایی همراه با آن بر گیاهان 
زراعی کاملا به شرایط محیطی بستگی دارد و ممکن است از محلی 
به محل دیگر تفاوت کند. از طرف دیگر مدل‌های گیاهی از جمله 
مدل ۷۷۵۳۵5۲ به‌دلیل دقت بالا در شبیه سازی ابزار مناسبی جهت 
پیش بینی عملکرد محصول در شرایط تغییر اقلیم است. هدف از این 
مطالعه بررسی اثر گرمایش جهانی از طریق کاربرد مدل گردش کلی 
۱ بر اساس سناریوهای انتشار 1360۳2.6 ,30۳45 و 
5 و مدل شبیه‌ساز رشد و ثمو ۱۷0۳05۴ در شرایط اقلیمی 
جنوب استان کرمان جهت شبیه سازی مراحل فنولوژی و عملکرد 
محصول سیب‌زمینی در شرایط پتانسیل می‌باشد. بنابراین» با توجه به 
اهمیت وضع هوا و اقلیم بر بخش کشاورزی» دستیابی به 
پیش‌بینی‌های بلند مدت اقلیمی در این بخش ضروری است 
(2019 ,ذه221). پیدایش مدل‌های اقلیم جهانی (0۷۳) و 
منطقه‌ای (50۷۲) این امکان ر به‌وجود آورده است تامحققان 
بتوانند تأثیر وضع هوا و اقلیم را در مقیاس‌های زمانی بلند مدت بر 
محصولات مختلف کشاورزی بررسی نمایند ( بلة 6۶ تات۳2ع۳۲ 
۵211 


با توجه به اهمیت گیاه سیب‌زمینی در جیره غذایی بشر و تاثیر 
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بیاتانی و همکاران. اثرات تغییر اقلیم بر عملکرد و مراحل فنولوژیکی گیاه سیب‌زمینی... ۶۷ 


پذیری آن از شرایط اقلیمی و محیطی بخصوص دما و رطوب ‏ در این 
تحقیق به بررسی آثرات تغییر اقلیم بر میزان عملکرد و طول دوره‌های 
رشد گیاه سیب زمینی پرداخته شد. به همین منظور از مدل 
ریزمقیاس‌نمایی آماری 57951 جهت ریزمقاس نمایی برونداد مدل 
گردش کلی جو ۵171182 تحت سناریوی های 1360۳ و مدل رشد 
و نمو ۷۷۵۳۵5۲ استفاده شد. نتایج حاصل از این تحقیق می تواند 
جهت اتخاذ روش های سازگاری با تغییر اقلیم مورد استفاده مسئولان 
سازمان جهاد کشاوررزی قرار گیرد. 


مواد و روش‌ها 

منطقه مورد مطالعه در این تحقیق ایستگاه‌های با طول دوره 
آماری ۷ ساله (۱۹۸۹ تا ۲۰۰۵) شامل جیرفت» کهنوج و منوجان می 
باشد. در این مطالعه جهت ریزمقیاس نمایی از داده‌های دمای حداکثره 
دمای حداقل» بارش, میانگین رطوبت نسبیء میانگین سرعت باد و 
ساعات آفتابی در مقیاس زمانی روزانه استفاده گردید. خروجی مدل- 
های گردش عمومی جو بزرگ مقیاس می‌باشند. بنابراین جهت ریز 
مقیاس نمایی خروجی مدل گردش عمومی از مدل ریزمقیاس‌کننده 
1 استفاده شد. در مدل 81957۷1 به‌منظور ریز مقیاس نمایی 
داده‌های پارامترهای اقلیمی از سه نوع داده الف- متغیرهای بزرگ 
مقیاس در دوره آینده (۲۰۱۱-۲۱۰۰) مدل گردش کی 
1 ب- متفیرهای بزرگ مقیاس ۱0 در دوره 
(۱۹۸۹-۲۰۰۵) و ج- داده‌های دیده‌بانی (دمای حداکثر. دمای حداقل» 
بارش, میانگین رطوبت نسبی, میانگین سرعت باد و ساعات آفتابی) 
روزانه مورد استفاده قرار گرفته است (2019 ,221201). به‌منظور 
تهیه داده‌های مزرعه‌ای جهت واسنحی و ارزیابی مدل ۷۷۵۳۵5۲ 
در شرایط پتانسیل (شرایط بدون محدودیت آب و غذاء کنترل بیماری- 
هاو علف‌های هرز) از داده‌های طرح‌های تحقیقاتی مرکز تحقیقات 
کشاورزی جیرفت طی سال‌های ۱۳۸۹-۹۳ جمع‌آوری و به‌عنوان دوره 
پایه یا دیده‌بانی مورد استفاده قرار گرفت. این اطلاعات شامل تاریخ 
کاشت» مراحل فنولوژیکی از کاشت تا سبز شدن» سبز شدن تا گلاهی 
و سبز شدن تا رسیدگی فیزیولوژیکی (برداشت) گیاه سیب‌زمینی در 
مناطق مورد مطالعه می‌باشد. به‌منظور کالیبراسیون مدل ۱۷۷۵۳۵5۲ 
با استفاده از داده‌های مزرعه‌ای از داده‌های فنولوژیکی و عملکرد در 
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سال‌های ۰۱۳۸۹ ۱۳۹۰ و ۱۳۹۲ و از داده‌های فنولوژیکی و عملکرد 
در سال‌های ۱۳۹۱ و ۱۳۹۳ برای تعیین اعتبار مدل استفاده شد. 


مدل تولید داده 515 

یکی از مهم‌ترین مدل‌های ریزمقیاس نمایی آماری مدل !5۳5 
است که مولد مصنوعی داده‌های هواشناسی می‌باشد. در این تحقیق. 
از داده‌های خروجی مدل گردش عمومی جو (20۳۴5[۷]2)) تحت سه 
سناریوی 0۳2.6 30۳4.5 و 13۳8.5 استفاده شده است. مدل 

41 در سال ۲۰۰۲ در انگلستان تهیه شده است. مدل 5۳0531 

روابط آماری را بر مبنای روش رگرسیون خطی چندگانه بین 
متفیرهای اقلیمی بزرگ مقیاس (پیش‌بینی‌کندهها) و محلی (پیش- 
بینی شونده‌ها) محاسبه می‌کند (2002 ,.21 66 ۷۷11107). در این مدل» 
قلیم محلی توسط اقلیم بزرگ مقیاس منطقه در فرم 1-۳0 بیان 
می‌شود. که در اینجا 1۶ نشان‌دهنده متغیر محلی است که کوچک 
مقیاس شده است» 2 محموعه‌ای از متغیرهای اقلیمی بزرگ مقیاس 
است و یک تابع تعیین مشروط به > است که بر اساس آموزش و 
اعتبارسنجی داده‌های تاریخی بدست می‌آید (2019 ,لصمات(هظ). 

در این تحقیق در ابتدا به‌منظور واسنجی و صحت یابی مدل 
متغیرهای بزرگ مقیاس مشاهداتی (جدول ۱) منطقه مورد 
مطالعاتی در مقیاس روزانه در دوره ۱۹۸۹-۲۰۰۵ (متغیرهای استخراج 
شده از ۲۳ ) به‌عنوان متغیر مستقل» و داده‌های دمای حداکثر» 
دمای. حداقل؛ بازش؛ میانگین رطویت نسبی: میانگین سرعت با و 
ساعات آفتابی روزانه مشاهداتی منطقه درهمین دوره. به‌عنوان متغیر 
وابسته به‌طور جداگانه وارد مدل شدند. پس از واسنجی مدل 81251۷۲ 
در ادامه به‌منظور بررسی توانایی مدل گردش کلی20۳:5[۷]2) در 
شبیه‌سازی متغیرهای اقلیمی منطقه, متغیرهای بزرگ مقیاس مدل 
2 در دوره ۱۹۸۱-۲۰۰۵ به مدل ٩1261۷]‏ (واسنجی شده 
در مرحله قبل) وارد گردیده و متغیرهای دمای حداکثر دمای حداقل» 
بارش, میانگین رطوبت نسبی, ميانگین سرعت باد و ساعات آفتابی 
منطقه برای این دوره کوچک مقیاس شدند. در نهایت. متغفیرهای 
کوچک مقیاس شده منطقه با متفیرهای مشاهداتی مورد مقایسه قرار 
گرفته و پس از کسب اطمینان از تونایی مدل 61551112 درکوچک 
مقیاس کردن متغیرهای اقلیمی منطقه. سری زمانی متغیرها با معرفی 
متغیرهای بزرگ مقیاس مدل 020۳1۷52 در دوره‌های (۲۰۴۰- 
۱ ۲۰۴۱-۲۰۷۰ ۲۰۷۱-۲۱۰۰) به مدل ۹151۷1 (واسنجی شده 


از مرحله قبل) برای دوره‌های آتی شبیه‌سازی شدند. 


۸ نشریه بوم شناسی کشاورزی. جلد ۱۶. شماره!. بهار ۱۴۰۳ 


جدول ۱- فهرست متغیرهای ۱۲۳ 
۷2۱۵۵6۵ ۱۲۲۷۳ 01 )115 -1 12۳016 


توضیحات 
جمناوتمعع(1 
سرعت باد زمینگرد در سطح ۵۰۰ هکتوپاسکال 
2 500 26 ۷۵۱0 «ملتتد منطومتادمع0) 
مولفه مداری سرعت باد در سطح ۵۰۰ هکتوپاسکال 
2 500 21 ۷/10 ۲۱01120۵121 
مولفه نصف النهاری سرعت باد در سطح ۵۰۰ هکتوپاسکال 
2 500 26 ۳10 70021 
اویی در سطح ۵۰۰ هکتوپاسکال 
2 500 21 ۲۷۵0۲007 
سمت باد در سطح ۵۰۰ هکتوپاسکال 
2 500 26 طمنامعن ۱۷۷180 
واگرایی باد در سطح ۵۰۰ هکتوپاسکال 
2 500 21 11۷۵۲۵۵86۵[ 
ارتفاع ژئوپتانسیل در سطیح ۵۰۰ هکتوپاسکال 
2 500 21 )طعنعط لمتامما000ع0) 
رطوبت نسبی در سطح ۵۰۰ هکتوپاسکال 
2 500 ]2 وانک‌تصط 1۵12۷۵ 
سرعت باد زمینگرد 
٩0۲]266 260500۱ ۷‏ 
مولفه مداری سرعت باد 
۵ )0۲120۵ همه 
مولفه نصف النهاری سرعت باد 
0 70021 ۱]۵66 
تاوایی 
۱ 
سمت باد 
ارت ۱۱۱۱۱۰۱۱ 
واگرایی باد 
٩۲]۲۵۵۵ ۵‏ 
سرعت باد زمینگرد در سطح ۵۰۰ هکتوپاسکال 
2 850 26 رهم۷۲۵۱ «ملتتد منطومتادمع0) 
مولفه مداری سرعت باد در سطح ۵۰۰ هکتوپاسکال 
2 850 26 ۷/180 ۲۱0۲20۵121 
مولفه نصف النهاری سرعت باد در سطلح ۵۰۰ هکتوپاسکال 
2 850 26 م1 70021 
اویی در سطح ۸۵۰ هکتوپاسکال 
2 850 26 ۲۷۵0۲67 
سمت باد در سطح ۸۵۰ هکتوپاسکال 
2 850 26 طمنامعتن ۱۷۷10 
واگرایی بد در سطلح ۸۵۰ هکتوپاسکال 
2 850 24 121۷۵۲۵۵86۵[ 
ارتفاع ژئوپتانسیل در سطلح ۸۵۰ هکتوپاسکال 
2 850 26 )طعنعط لمتاجما000ع0 
رطوبت نسبی در سطح ۸۵۰ هکتوپاسکال 
2 850 26 بوانک‌تصط 1۵12۷۵ 
میانگین فشار سطح دریا 
اوعد 1۵۷۵ 562 م۷۲62[ 
رطوبت نسبی نزدیک سطح زمین 
تیه وتاقاه؟ معویه تععل 
رطوبت ویژه نزدیک سطح زمین 
تصیط مرللمعمو معوناه بقل 
دمای هوا در ارتفاع ۲ متری 
هرهاق 2-0 


۱ 


ره 
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5-7 


۳3 7 


۱ 


۳82 


۳950 


1850 
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ردیف 
10۳ 
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بیاتانی و همکاران. اثرات تغییر اقلیم بر عملکرد و مراحل فنولوژیکی گیاه سیب‌زمینی... ۶۹ 


معرفی مدل ۲۷۷۵۲۵5 

مدل 170۳057 نمو فنولوژیکی و رشد گیاه را از آغاز جوانه زنی 
تا زمان رسیدگی بر اساس خصوصیات ژنیتیکی گیاه و شرایط محیطی 
شاشن کی شون یه ترا تادی کک سرعت تفای حلت 
02 توسط کانوپی استء که به انرژی تابشی جذب شده بستگی 
دارد که تابعی از تشعشع روزانهه سطح برگ گیاه و ضریب خاموشی 
نور در کانوپی است. ۱۷۵۳05 یزان جذب 002 اخالص از ینک 
گیاه را بر پایه تششع فصال فتوسنتزی (۵۰ درصد از کل تشعشع 
خورشیدی با طول موج ۴۰۰ تا ۷۰۰ نانومتر) جذب شده توسط گیاه 


83 3 2 (ونامعع۱۵۲ اصاونا 


شبیه سازی می‌کند. 002 با استفاده از انرژی تأمین شده بوسیله نور 
به کربوهیدرات (۳120) تبدیل می‌شود. قسمتی از کربوهیدرات‌های 
باقیمانده در مسیر تنفس رشد به ماده خشک تبدیل می‌شوند. 
فرآیندهای ذکر شده ساختار مدل در شکل ۱ ارائه شده است. رشد در 
مدل به‌صورت دینامیک با گام زمانی یک روزه شبیه سازی می‌شود ( 
2 باه به تصمتنصع؟ 122 مسدل ۷۷۵۳۵۹۲ دارای 
نسخه‌های مختلفی است که در طول سال‌های مختلف ایجاد شده 
است. در اين مطالعه از نسخه 7.1 ۷۵۲ ۷۷۵۳۵۹۲ استفاده شد. که 
توسط مرکز تحقیقات اگروبیولوژی و حاصلخیزی خاک و مرکز 
تحقیقات کشاورزی هلند در دانشگاه واخنینگن ساخته شده است. 
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شکل ۱- مراحل شبیه سازی رشد و نمو گیاه در مدل ۲۷۵۳۵5۲ (1994 وله اه اذمک) 
(1994 ولو )6 )6۵) ۵۵06 ]۷۷۵۲۵۹ 1 موه 20۲۵ کصعان 01 دععد)5 -1 ۲۱۵۰ 


برای توصیف نمو فنولوژیکی گیاه در مدل 9/0057 از متفیر 
5 استفاده می‌شود که برای جوانه زنی» گلدهی و برداشت به 
ترتیب مقلدیر ۰» ۱ و ۲ در نظر می‌گیرد که تابعی از درجه روز می‌باشد 
(2014 ,لد اه نصه‌تنصهگ 126). دوره زندگی سیب‌زمینی در مدل 
به دو مرحله تقسیم می‌شود: 

جوانه زنی (0۷5-0) تا زمان گلدهی (۷5-1) 

دوره زایشی: از گلدهی (12۷۹-1) تازمان رسیدگی 
فیزیولوژیکی (0۷5-2) 


برای اجرای مدل ۱/۸ الف) به داده‌های هواشناسی 
متوسط سرعت باد و بارش) به‌صورت روزانه و ب) داده‌های مزرعه‌ای 
لازم می‌باشد. داده‌های مورد نیاز برای کالیبراسیون و تعیین اعتبار 
مدل از آزمایشات مزرعه‌ای انجام شده در مرکز تحقیقات کشاورزی 
جیرفت فراهم شد. مرحله کالبیراسیون مدل ۷۷۵۳05۲ با استفاده از 
متغیرهای اندازه‌گیری شده رشد گیاه زراعی نظیر مشاهدات مراحل 
فنولوژیکی گیاه سیب‌زمینی (کاشت تا سبز شدن» سبزشدن تا گلدهی 


۰. نشریه بوم شناسی کشاورزی. جلد ۶ شماره!. بهار ۱۴۰۳ 


و سبزشدن تا رسیدگی فیزیولوژیکی (برداشت)) و میزان عملکرد در 
سال‌های ۰۱۳۸۹ ۱۳۹۰ و۱۳۹۲ انجام گرفت. بعد از کالبیراسیون 
نمودن مدل ۰۷۷۵۳0۵5 تعیین اعتبار مدل با استفاده از داده‌های 
مزرعه‌ای سال‌های ۱۳۹۱ و ۱۲۹۳ انحام گرفت. 


معیارهای ارزیابی عملکرد مدل 

جهت ارزیابی توانمندی و دقت مدل ٩۲25[۷‏ در شبیه‌سازی داده- 
های دیده‌بانی (۱۹۸۹-۲۰۰۵) و همینطور تعیین اعتبار و ارزیابی مدل 
۲ در شبیه‌سازی مراحل فنولوژیکی مختلف گیاه سیب- 
زمینی و عملکرد غده سیب‌زمینی از معیارهای آماری» همچون ریشه 
دوم میانگین مربعات خطا (1۷]5۲)» جذر میانگین مجذور خطای 
نرمال شده (151907)» حداکثر خطا (۷10» ضریب تبیین ( 
ضریب ناش ساتکلیف ( :م۳) استفاده گردید (2019 ,نعم/ه۵). 


۳( ]1 - )یج 2116۴8 


معادله (۱) 

معادله (۲) .| ,۲۸ - ,۲۵25 < ۷18 
۱ 

معادله (۳) مور - رورم 


در این معادلات ٩‏ و ]۷ به‌ترتیب مقادیر شبیه‌سازی و اندازه- 
گیری شده. 1 میانگین مقادیر 1۷1 و 2 تعداد مشاهدات می‌باشد. 
مقادیر جذر میانگین مجذور خطا و جذر میانگین مجذور خطای نرمال 
شده در حالت اپتیمم یا حالتی که مقادیر شبیه سازی شده و 
اندازه‌گیری شده مساوی باشد برابر با صفر و مقدار کارایی مدل در 
ای لیوا یکتم ناف 

همچنین می‌توان برای ارزیابی مدل ۷۷۵۳05۲ از نمودار مقادیر 
واقعی (۷11) در برابر مقادیر شبیه‌سازی شده (51) و از معادله رگرسیون 
خط ۱:۱ رسم شده استفاده کرد. 
معادلعه (۴) (85,)5 - ,6۲ ع الا 
۰ کم ۳ 


سپس اختلاف عرض از مبدا (؟ ۳ 


همان شیب خط رگرسیونی به‌ترتیب با صفر و یک مورد مقایسه قرار 


که 


۶۲ 0۲۲۵ 20ع]۱۷ ]۴۵0 -1 
-۲۲۲۵۲ عتهیان5 1۷۲22 00۲ -2 
۲ ۱۷۲231۳۲۷۱ -3 

4 - 15-۹00۵۲60 ۲ 

تلهم 11116ماناه وهآ - 5 


کالییر اسیون و اعتبار سنجی مدل رشد گیاهی ۲۷۷۲۵۲۵۵۲ 
تحت شرایط مزرعه ای 

پارامترهای مختلف مدل رشد گیاهی ۷۷۵۳0۵5۲ که با هدف 
حداقل کردن اختلاف بین مقادیر اندازه‌گیری شده و شبیه‌سازی شده 
عملکرد غده و طول مراحل فنولوژیک سیب‌زمینی در مرحله 
کالیبراسیون نهایی شد در جدول ۲ آمده است. از بین این پارامترها 
1 و ۲51/۸2 دمای پایه رشد 151110511 سطح ویژه 
برگ, حداکثر سرعت جذب 702 بهره وری تبدیل به برگ» ريشه و 
ساقه طبق جدول ۲ تعدیل شدند و سایر پارامترها به‌صورت مقلایر 
پیش فرض مدل در نظر گرفته شدند. 

بر اتباتی کایج سقاذیر انتازهگری شیه رت #سال بتارم 
کاشت‌های مرسوم (جدول ۲) میانگین طول دوره کاشت تا سبز شدن 
۷ روز ميانگین دوره سبز شدن تا آغاز گلاهی ۳۴ روز و میانگین 
طول دوره سبز شدن تا رسیدگی فیزیولوژیکی (برداشت) ۱۰۲/۷ روز 
بطور متوسط برای منطقه مورد مطالعه محاسبه و تعیین گردید. مدل 
در مرحله کالیبراسیون طول دوره سبز شدن را ۲۴/۳ روز دوره آغاز 
گلاهی را ۲۳/۷ روز و طول دوره سبز شدن تا رسیدگی فیزیولوژیکی 
را (برداشت) را ۱۰۰روزء بطور میانگین برای ۳ ایستگاه مورد مطالعه 
شبیه‌سازی نمود. در مرحله تعیین اعتبار مدل بطور میانگین طول دوره 
سبز شدن را ۲۳/۷ روز طول دوره گلدهی را ۳۴/۲ روز و طول دوره 
سبز شدن تا رسیدگی فیزیولوژیکی (برداشت) را ۱۰۱/۸ روز برای 
تام گنای مورخطاله سای نیو 

مقایسه عملکرد شبیه‌سازی شده با مقادیر واقعی در مرحله 
کالبیراسیون و اعتبارسنحی در شکل ۲ نشان داده شده است. با توجه 
به شکل ملاحظه می‌شود که مقدار ضریب 1 محاسبه شده برای 
مراحل فنولوژیک و عملکرد در هر دو مرحله کالیبراسیون و اعتبار- 
سنجی به‌ترتیب بیش از ۰/۸۷ و ۰/۹۴ بود که نشان‌دهنده دقت بالای 
مدل در شبیه‌سازی مراحل فنولوژیک و عملکرد سیب‌زمینی است. 

نتایج تعیین اعتبار مدل ۷۷0۳51 برای مراحل فنولوژیک و 
عملکرد سیب‌زمینی در جدول ۴ نشان داده شده است. مدل 
۲ برآورد قابل قبولی برای هر دو متغیر به همراه داشت. بر 
اساس نتایج حاصل از شبیه‌سازی مشخص شد میزان عملکرد شبیه- 
سازی شده. برازش خوبی با مقادیر دیده‌بانی شده دارد (شکل ۲). 


بیاتانی و همکاران. اثرات تغییر اقلیم بر عملکرد و مراحل فنولوژیکی گیاه سیب‌زمینی... ۷۱ 


جدول ۲- پارامترهای گیاهی مدل ۲۷۵۲5۲ در مرحله کالیبراسیون در شرایط یتانسیل 
مصممصم تامهم 06۲ص مها ممته ای ما ۵۵08 ۷۷۵۲۵۵ ۲66 ۵ وامام2 2۲ ۳۱۵۸۲ -2 12016 


مقدار 
واحد ۲ 
انا ۱/۹ توصیف پارامتر پارامتر 
تعدیل شده پیش فرض جمنارننعع0 ۵9۳۵۵)۵۲ ۳2۲ ۲( 
اششاااان ۱۱۱ تاالت ۱۷۸۵ 
۹ 8 10۲ ۲۵5۲010ظ) مصا ۵0100 وعتتتاج تمصع[ 
درجه روز ودوه 170 و روز از کاشت تا سبز شدن ۱ 903/0۳2 
۹4 6 0۵ ع )مهن 1۳0 5۵1 068760 انعر[ 
درجه روز 501 150 ۱ مجموعه درجه روز از سبز تن گلاهی ۱ ۳۹0/1 
4« 1107 م) معصوعتمصصی حطم ممتامع11مع مصصت) وح(1[ 
یت مجموعه درجه روز از گلاهی تا بلوغ یا برداشت 
271 2 59154 1530 ۶۲ اتمه ما عصته 1107 1۳0 اعو 068166 122717[ 1-32( 
۵ 
هکتاراکیلوگرم 0002 0003 ۳ ٩۸8‏ 
٩060121 ٩۱1۴2606 0۴ 5 21"‏ 
گرم /هکتار /اساعت اکن مه 
35 30 ت ای بمع ون ۱۹ 
اسبط مدای علماود دمن 0 عاق حصتاصدن‌نه]۷ 
5 5 1 گ 
کیلوگرم اکیلوگرم 00 0120 ۰ بهرهوری تبدیل به بر "0 
ص1۶2 ۷ 60۳0۷۵۳5109۱ 62[ 
کیلوگرم/کیلو گرم و055 0220 بهره‌وری تبدیل به ريشه 
.ی ۲ 36 ۵00۴65 اد متا 0ترظ 


جدول ۲- طول دوره مراحل فنولوژیکی بر اساس تاریخ کاشت (تعداد روز از اول زانویه) بصورت میانگین طی سال‌های ۱۳۸۹-۹۳ 
2011-15 ۷6۵۵۲5 6ظ) ال هه مه (1 متا صول صمظ بقم) 026۵ مصصهام چم 0ععدوها مها ام‌نعماممهمطم ۵1 «متاه۳۳ -3 12016 
مراحل فنولوژیکی (روز) 
ک6ع۶]2 اهع1ع۳۵۵۸۵۱0 
سبز شدن تا رسیدگی 


فیزیولوژیکی (برداشت) تسیز سدن ۷ کوج تاریخ کاشت شهر 
آمنعمامنووطم ت50 ۱ 0 سال دن‌ 
1 عصاصهام ]0 206 
۲ 9۱۵ (02) (ومة) 9۳( 
دیدبانی شبیه سازی دیدبانی شبیه سازی دیدبانی شبیه سازی 
جمتاججعهعوای .۰ 0عاقاصت... حمتاعهوهام ...۰ 0عاقات.. ممتاهعهووام ...۰ 60اقلنامهزه 
106 104 37 35 26 25 20 1389 
100 98 34 33 24 23 20 1390 مخ 
102 100 37 36 26 24 20 1391 و 
105 103 36 34 25 24 20 1392 پانال 
10 106 35 33 24 22 20 1393 
103 104 33 32 24 24 20 1389 
101 102 36 35 25 24 20 1390 ۲ 
103 104 36 34 24 22 260 1391 ۱ کهنوج 
98 96 33 32 23 22 20 1392 [هتاعطفگ 
102 101 35 34 26 24 260 1393 
104 102 35 33 22 21 20 1389 
99 91 35 34 24 23 20 1390 دا 
103 103 35 33 22 20 20 1391 ات 
101 100 33 32 22 21 260 1392 مدزمنمه]۷ 
105 104 36 35 23 22 20 1393 
10۳2.7 1016 35 2317 24 227 ست ست ميانگین 


(0 


۲ نشریه بوم شناسی کشاورزی. جلد ۱۶. شماره!. بهار ۱۴۰۳ 


مقادیر ميانگین مربعات خطا برای هر دو متغیر در محدوده عالی 
قرار داشت 1 از ۱/۸ درصد تجاوز نکرد» بعالاوه شیب خط رگرسیون 
بین مقادیر دیده‌بانی و شبیه‌سازی شده برای مراحل فنولوژیک و 
عملکرد با شیب خط ۱:۱ تفاوت معنی‌داری نشان نداد (زشکل۲). بر 
اساس شکل ۲ بیشتر از ۸۷ درصد مراحل فنولوژیک و بالای ٩۰‏ 
مراحل فنولوژیک و عملکرد دیده‌بانی شده قرار دارد. همچنین ضریب 


تبیین () نیز همبستگی بالایی (0.87< *) بين مقادیر دیده‌بانی 
شده با مقادیر شبیه سازی شده نشان می‌دهد. مقادیر 1۷1۳ نشان‌دهنده 
(12 >۷]۳) کارایی بالای مدل می‌باشد. بنابراین» نتایج بررسی و 
مقایسه شاخص‌های مختلف آماری نشان‌دهنده دقت بالای مدل رشد 
گیاهی ۷۷0۳057 در شبیه سازی عملکرد و مراحل فنولوژیک 


سیب‌زمینی است. 


جدول 4- نتایج کالیبراسیون و تعیین اعتبار مدل برای عملکرد و مراحل فنولوژیک سیب‌زمینی در شرایط پتانسیل 
وصمانصم اهناصماهمه 0ص مع‌هه)ه ماع ه۵امصعمطامر مه ۱۵۱0 2)0)هه رم مت)۷۱۵۱۱02 061مصظ مصه وعاوع ماه نان -4 م1201 


۷ نع طف۳ 1۳/98 بر 
4 093 4 1 24 
ور 4و0 2 05 23.5 
44 1 4 09 218 
7 075 49 137 35 
58 081 4 13 33 
3 0.15 2.5 02 33 
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5 ۰ 072 ۳7 17 100 
7 095 17 1۳ 997 
8 ۰ 097 25 5675 445 
5 099 2.5 5675 41 
2.6 02 1.6 32906 ۹0۹4+( 


دیده‌بانی ابستگاه مراحل فنولوژیک و عملکرد 
٩) 0‏ جمناهاگد ‏ عععهاه افتعمامممصام 0و ۱۷۱۵۱0 
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بیاتانی و همکاران. اثرات تغییر اقلیم بر عملکرد و مراحل فنولوژیکی گیاه سیب‌زمینی... ‏ ۷۳ 
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شکل ۲- مقایسه مقادیر شبیه‌سازی شده با مقادیر اندازه‌گیری ده برای عملکرد غده گیاه سیب‌زمینی (داده‌های سال‌های ۰۱۳۸۹ ۱۳۹۰ و ۱۳۹۲ 
برای مرحله‌ی کالبیراسیون و داده‌های سال‌های ۱۳۹۱ و ۱۳۹۳ برای مرحله‌ی تعبین اعتبار)؛ خط یک به یک نیز در شکل نشسان داده ده است. 
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٩۳0۱۷۱ 1۱ ۵ 6‏ 0و2 وز مصنا عصم-ه)-عصه معط ورصم)02ناه ۷ هد 6ظ) 10 2015 280 2013 10۲ 0262 0ص2 ععهاو صمتاه تنل 


تولید سناریوهای تغییراقلیم در دوره‌های آتی 

عملکرد مدل ٩1251۷‏ در ریزمقیاس‌نمایی داده‌های اقلیمی 
وسیله‌ی شاخص‌های آماری (جدول ۵) مورد ارزیایی قرار گرشت, 
قابل توجهی برخوردار است (بالای ۰/۹۹). حداکثر خطا (015) برای 
شبیه سازی دای حداقل» دمای خداکر و بارش بظور مانگین رای 
هر سه ایستگاه مورد مطالعه به‌ترتیب ۰/۱۷ ۰/۱۵ و ۱۹/۹ بدست آمد 
که نشان از خطای بیشتر شبیه سازی بارش نسبت به پارامتر دما می 
باشد. شاخص ریشه میانگین مربعات خطا (1۷151) در شبیه سازی 


بارش نیز نسبت به دیگر پارامترها به‌مراتب بیشتر بوده (۲/۲) و خطای 
پیشتری را نشان می‌دهد. بنابراین, با توجه به نتایج جدول ۵ توانایی 
مدل ٩۲95۷‏ در شبیه سازی پارامترهای هواشناسی تایید شده. ولی 
با این حال مدل در شبیه سازی بارش دقت کمتری را نشان می‌دهد 
که با یافته‌های جعفرپور 2016 ,7ا1212700)( نظری و همکاران 
(21.,2015 6 22271) و طایی سمیرمی و همکاران 12:6) 
(2014 ,.له اه مصصم‌تنصدگ مطابقت دارد. در ادامه نتایج بدست آمده 


از ارزیابی مدل 51751۷ در جدول ۵ ارائه شده است. 


بیاتانی و همکاران. اثرات تغییر اقلیم بر عملکرد و مراحل فنولوژیکی گیاه سیب‌زمینی... ‏ ۷۵ 


جدول ۵- کارآیی مدل 51051 در شبیه‌سازی بارش و میانگین دمای حداقل و دمای حداکثر در دوره دیدبانی با داده‌های ۷0۲۴ 
0 ۵۵۵0۵60۱۲۵ حصا یه 20 صمرتصه همه له صمااه)نمرتع 0۳ عصت/مامسته دز ع0مصه 6 ع.) ۵۶ وم‌جهنم ۳ -5 12016 
2 ۱۱۲ ۱۱۲۵ ۵۵0 صمناه ۲ ووهاه م6 مس 68۵۵۵۵6 ۱۸22311001۵ 


معیار دقت 
متماتن و۸ 


1۷/9 ۰ (160 )60( 11۳ 82 
002 01 021 1 
003 01 024 1 

وو0 21 0 12 
1 02 01 002 
1 01 005 002 

وو0 108 04 07 
1 01 001 001 
1 01 01 003 

وو0 28 08 16 


در جدول ۶ نتایج مربوط به شبیه‌سازی دمای حداقل» حداکثر و 
بارش در سه دوره آینده تحت سناریوی انتشار 30۳2.6 ,0۳4.5 و 
5 منمایش داده شده است. با توجه به جدول ۶ دمای حداکثر و 
حداقل شبیه‌سازی شده برای دوره آینده (۲۰۱۱-۲۱۰۰) تحت 
سناریوهای 1۳ نسبت به دوره دیدبانی در هر سه سناریو 50۳ 
افزايش نشان می‌دهد. بارش نیز در دوره‌های آینده در هر سه 
سناریوها دچار تغییر شده است و روند روبه کاهش را نشان می‌دهد. 
میانگین افزايش سالانه دمای حداقل و حداکثر در مناطق مورد مطالعه 


متغیر ایستگاه 
د, صد خطا دا ۱/2 0 
۳ ِ 
۴ ۳۵۲۵۵۵۵۵۵6 
و0 02- دمای خبافل 
طلصص ]1 
مای حداکة 
099 002 ای ثِ#- 
-[ 1 4 
و0 45 بارش 
ارفافاتهازه ات ۱9 
و0 003 دمای خباال 
صتص- ]1 
5 7 
و0 0004 وا م 2 
-[ 1 2 
9 
- 
وو0 7 بارش 
ارفافاتهازه ات ۱9 
099 04 دماق خاي 
صتص- ]1 
مای حداک: کِ 
و0 004 ای ام 
1-2 4 5 
۱ 
[ 
099 1.8- ۱ بارش 
ارفافاتهازه ات ۱9 


از ۰/۴ درجه سانتیگراد در سناریوی 10۳2.6 تا ۰/۷ درجه سانتیگراد 
در سناریوی 130۳8.5 مشاهده شد. بیشترین افزایش دمای حداکثر و 
حداقل به‌ترتیب به‌میزان ۱ و ۱/۳ درجه سانتیگراد درایستگاه کهنوج و 
جیرفت در سناریوی 5۳8.5 مشاهده گردیدد. میانگین کاهش بارش 
به‌میزان ۱/۲ میلیمتر در سناریوی 5۳2.6 تا ۱/۷ میلیمتر در 
سناریوی 10۳8.5 نیز مشاهده گردیدد. بیشترین کاهش بارش به 
مقدار ۳۸ میلیمتر در ایستگاه جیرفت مشاهده شد. 


۶ نشریه بوم شناسی کشاورزی. جلد ۶ شماره!. بهار ۱۴۰۳ 


جدول 7- میزان تغییرات میانگین دمای حداقل و حداکثر و بارش سالانه در ایستگاه‌های منتخب در سه دوره ۲۰۱۱-۲۰۰ ۲۰۱-۲۰۷۰ و 
۲۰۷۱-۶۰ نسبت به دوره دیدبانی تحت سه سناریوی انتشار 110۳2:6 10۳4.5 و 10۳8.5 


6 ۱ مصمناهاه ۹0۱60660 ۵۲ مماهازمنمع 0۳ احبصصه مضه ده مره تفه مضه منم 2۲۵2۵6 66 صز وم‌عصمطان -6 م1201 
۵۵۵۵05 009510۱ ۳۲۵۵ )6 06ص مه مناج عوهان عمط 6۵ 0ع۲دصمع 20100-2071 290 2070-2041 ,2040-2011 08من۵ع ۲6۵ 
۹5 /20۱ 1۸۳4.5 ,1۸0۳2۰6 


بارش دمای حداکثر دمای حداقل 
)۳۵01 6ج مرصه) صصتصنآ ۱ سرصه) صصتصنآ ۱ ایستگاه 
منوجان کهنوج جیرفت منوجان جیرفت منوجان کهنوج جیرفت «مناهلا 
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۱0۳-٩1 5 1 0.4 1 03 -8 -4 7 -4 9 -18‏ 2 0.4 2 04 2 1۹208 5 
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بررسی تأثیر تغییر اقلیم بر طول مراحل فنولوژیکی و 
عملکرد غده 

تلفیق نتایج حاصل از مدل اقلیمی 020۳5(۷/2) تحت سناریوهای 
انتشار 10۳2.6 ,30۳45 و 30۳8.5 در سه دوره آقلیمی ۲۰۴۰- 
۱ ۲۰۴۱-۲۰۷۰ و ۲۰۷۱-۲۱۰۰ با مدل شبیه‌ساز رشد و نمو 
7 تشان داد که طول مراحل مختلف فنولوژیکی, از کاشت 
تا سبزشدن» سبزشدن تا گلاهی و سبزشدن تا رسیدگی فیزیولوژیکی 
(برداشت) (بدون محدودیت آب و عناصر غذایی) و عملکرد غده 
سیب‌زمینی در شرایط پتانسیل در مناطق مورد مطالعه کاهش خواهد 
یافت. نتایج حاصل شده با یافته‌های ادواری و همکاران ( 66 ۸0۷271 
90 ,.21) و ترکمان (2016 ,حقص1:1) مطابقت دارد. 

بر اساس نتایج شبیه‌سازی در مناطق مورد مطالعه در دوره 
اقلیمی آینده. طول دوره کاشت تا سبز شدن بطور میانگین بین ۰/۳ 
روز تا ۱۸۲ روز کاهش خواهد یافت. در مناطق گرمسیری که دارای 
زمستان‌های معتدل هستند و کشت پاییزه سیب‌زمینی انجام می‌شود. 
مرحله‌ی کاشت تا سبز شدن بیشتر از سایر مراحل فنولوژیکی سیب- 
زمینی. تحت تأثیر اثرات ناشی از گرمایش جهانی خواهد بود. البته تا 


کنون تأثیر منفی کاهش طول دوره کاشت تا سپزشدن به‌علت 
گرمایش جهانی بر عملکرد نهایی غده‌های سیب‌زمینی گزارش نشده 
است (2016 ,27020 ]). با اینوجود افزایش دماء ناشی از 
تغيرقلیم طول دوره رشد را در مناطق گرمسیری به‌شدت کاهش 
خواهد داد. بنابراین» موضوع تغییر تاریخ کاشت یکی از آثار گرمایش 
جهانی و راهکار سازگاری با پدیده گرمایش جهانی می‌باشد 
(2016 ,صعصصیت). طول دوره سبزشدن تا گلدهی بطور میانگین 
بین ۰/۷ روز تا ۱/۵ روز کاهش خواهد یافت. بطور متوسط میزان 
کاهش دوره‌ی گلدهی در سناریوی ۳2.6 بین ۰/۵ تا ۰/۸ روز در 
سناریوی 50۳4.5 در محدوده‌ی ۱ تا ۱/۶ روز و در سناریوی 
5 در محدوده‌ی ۱/۰۵ تا ۱/۷ روز خواهد بود. مرحله‌ی 
سبزشدن تا رسیدگی فیزیولوژیکی غده‌ها (برداشت) به‌علت طولانی‌تر 
بودن آن بیشترین کاهش طول دوره در اين مرحله اتفاق فتاد. به‌طور 
متوسط میزان کاهش طول دیره رسیدگی فیزیولوژیکی غده‌ها 
(برداشت) برای ۳ ایستگاه مورد مطالعه بین ۱/۵ روز در سناریوی 
6 تا ۳/۶ روز در سناریوی 6۳8.5 قرار دارد (شکل ۳). 
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بیاتانی و همکاران. اثرات تغییر اقلیم بر عملکرد و مراحل فنولوژیکی گیاه سیب‌زمینی... 
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شکل ۲- شبیه‌سازی مدت زمان کاشت تا سبز شدن (بالا سبز شذن تا گلدهی (وسط) و سبز شدن تا رسیدگی فیزیولوژیکی (برداشت) (پایین) 
گیاه سیب‌زمینی بر اساس مدل 027۳:5۷۲2 برای سه دوره اقلیمی ۲۰۱۱-۳۰۰ ۲۰۱-۲۰۷۰ و ۲۰۱۷۱-۳۲۱۰۰ تحت سناریوهای 10۳ با 
استفاده از مدل ۲۷۵۲۵5۲ 


۲ (۱)2ه۱0م۱و دض م6 ععصمو مه مه انم مصته 110۲ م6 ماه رروه)) ععصععمصه ۲۵ مصصنا عصتاصدام ۵ ممتا بصن -3 ۲۱۵۰ 
۲ 2100-1۴ ۸۳0 2070-2041 ,2040-2011 ۵۲۱005 16اجصصنا عع۲ظ) 10۲ ۵061 02۵6۱۷۲۸2 ون 2960 اصهاص مهم 0۴ رصصم)00۲) 
۱۵۵ ۷۲۷۵۲۵۹۲ 66 هون وماتفصهعو 1۱0۳ 


شبیه‌سازی عملکرد غده سیب‌زمینی 

بر اساس گزارشات موجود واکنش عملکرد سیب‌زمینی نسبت به 
تغییيراقليم بسته به منطقه و سناریوی تغییر اقلیم متفاوت است 
(2016 م,ححصصیتن]). نتایج شبیه‌سازی عملکرد غده سیبزمینی 
نشان‌دهنده‌ی کاهش عملکرد در هر سه ایستگاه مورد مطالعه است 
(جدول ۷ با توجه به تتایح بدست آمده. بیشترین میزان کاهش 
عملکرد در دوره سوم (۲۰۷۱-۲۱۰۰) مشاهده شد. درصد کاهش 
عملکرد غده بطور میانگین بین ۲ تا ۴/۲ درصد در مناطق مورد 
مطالعه مشاهده شد(جدول ۷). 
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در ایستگاه جیرفت در سه دوره آتی بطور متوسط تحت سناریوی 
6 ۵ فرص ار یی 5 ۲/۱۰۵۳ حرف وس برس 
5 ۴/۰۵ درصد کاهش عملکرد مشاهده گردید. در ایستگاه 
کهنوج و منوجان به‌ترتیب در سناریوی 5۳2.6 ۱/۶ و ۲/۹ درصد. 
در سناریوی 50۳4.5 ۲/۹ و ۵/۲ درصد و در نهایت» سناریوی 
05 ۴ و ۴,۸ درصد کاهش عملکرد غده مشاهده شد. بنابراین» 
با توجه به نتایج بدست آمده بیشترین کاهش عملکرد غده سیب‌زمینی 
در هر سه دوره آینده تحت سناریوهای اقلیمی 3۳ در ایستگاه 


منوجان می‌باشد. 
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۸ نشریه بوم شناسی کشاورزی. جلد ۶ شماره!. بهار ۱۴۰۳ 


جدول ۷- شبیه‌سازی عملکرد غده سیب‌زمینی (کیلوگرم در هکتار) و درصد کاهش عملکرد بر اساس مدل 20۳/51۷2 تحت سناریوهای 130۳ در 
مناطق مورد مطالعه در دوره‌های ۲۰۱۱-۲۰2۰ ۲۰۱-۲۰۷۰ و ۲۰۷۱-۲۱۰۰ 
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اثرات تغییرات مدت زمان مراحل فنولوزیکی بر عملکرد 
سیب زمینی 

به‌منظور ارزیابی روند عملکرد سیب‌زمینی در مناطق مورد مطالعه 
تحت تأثیر تغییرات ایجاد شده در طول دوره رشد گیاه سیب‌زمینی 
(کاشت تا سبزشدن» سبزشدن تا گلدهی و سبزشدن تارسیدگی 
فیزیولوژیکی (برداشت)» نتایج حاصل از شبیه‌سازی دمای حداقل و 
دمای حداکثر با نتایج شبیه‌سازی مراحل فنولوژیکی و عملکرد غده 
سیب‌زمینی تلفیق گردید. نتایج بیانگر افزايش میانگین دماهای حداقل 
و حداکثر به‌میزان ۰/۳۷ تا ۰/٩‏ درجه سانتیگراد طی سه دوره آینده 
می‌باشد (جدول ۸). طول مراحل فنولوژیکی گیاه سیب‌زمینی از سبز 
شدن تا رسیدگی فیزیولوژیکی (برداشت) به‌موازات افزایش میانگین 
دمای حداقل و حداکثر کاهش یافته‌اند (جدول ۸). کاهش طول 
مرحله مییاشدن بین 9/۲ ۱/۹ رون رخا گلذهی‌بیین ۰/۵ ۱۸ 


روز و مرحله سبزشدن تا رسیدگی (برداشت) بین ۷/۴ تا ۳/۵ روز قرار 
دارد. نتایج جدول ۸ حاکی از آن است که با افزایش میانگین دما 
طول مراحل رشد فنولوژیکی گیاه کاهش خواهد یافت و به تبع آن 
کاهش میزان عملکرد غده سیب‌زمینی را در پی خواهد داشت. 
همانطور که در جدول ۸ نیز مشاهده می‌شود میزان کاهش عملکرد 
غده در سناریوی 5۳2.6 در دوره اول (۲۰۱۱-۲۰۴۰) از ۶۰۳/۶ 
کیلو در هکتار به‌مقدار ۲۵۹۲/۳ کیلو در هکتار در سناریوی 10۳8.5 
در دوره سوم (۲۰۷۱-۲۱۰۰) رسیده است. بیشترین کاهش در طول 
مراحل فنولوژیکی مربوط به مرحله سبزشدن تا رسیدگی (برداشت) به 
میزان ۳/۵ روز است. شکل ۴ نشان می‌دهد که به‌ازای هر یک درجه 
سانتیگراد افزايش میانگین دما سالانهه عملکرد غده به‌میزان ۴/۹۹ 
درصد کاهش خواهد یافت. 


جدول ۸- مقایسه مدت زمان کاهش مراحل فنولوژیکی (روز) و عملکرد غده سیب‌زمینی (کیلوگرم در هکتار) و افزایش میانگین دمای حداقل و 
حداکثر (درجه سانتیگراد) بر اساس مدل گردش عمومی 017۷]52) تحت سناریوهای 13۳ در مناطق مورد مطالعه 
6۵ 2 (2 / 16۵) ۷1۵۱0 مهاب ماهامض مضه (وودل) هه اهمزوماممهمصان ۵ دمتام ۳۵0 ۵۴ صمتاح سل عمط ۵۶ مموزبدرصمی -8 12010 


۱/۸ 261 مظ صمناهاین تمصع عمط چم ود (فستواعن) 06۵۲665) همه یمه 0ص صیصصتص ۸۷۵۲۵۵۵ ۵1 
۵۵ ۲007 6ظ۳) سا دمدفصمعه ۱ ع۵صت 


2071-10 2041-00 2011-40 پارامترها 
5 ۳4,5 10۳2.6 ,10۳8 ۳0۳4.5 . 10۳2.6 0۳8,56 . 10۳4.5 1۸0۳2,6 ۱ 
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۰.3 2077.8- 94 4 ۰ ۰ 1( 6.3 21- 38- 036- 0 
09 07 04 08 0.6 005 038 037 037 0۳۵ ۸۷۵۲۵۵6 


بیاتانی و همکاران. اثرات تغییر اقلیم بر عملکرد و مراحل فنولوژیکی گیاه سیب‌زمینی... ۷۹٩‏ 


1 1.5 


۲7 < 4.99617 7 1 


18*2 -<- 6 


(6۵۲۵2560 ۵۷۵۵۲ ۷۱۵۱0 )۷۸( 


هه ۸۲۵۲۵۵ 


شکل - درصد کاهش عملکرد غده سیب‌زمینی به‌ازای افزایش میانگین دمای سالانه با استفاده از مدل گردش عمومی 020۳۷۲52 تحت 
سنار یوهای 1560۳ در مناطق مورد مطالعه 
امومع ۵۵۳۱۷۲۹2 فصو ده مره امسصصه مفو ماه صد مومع مها رم ۷۱۵۱0 وهای ماهامص صز »وه6 068۵ ۳۵۲۵۵۵6۵۸۵6 -4 .۲12 
۵ /۹0۷۱6 ۲۵۵ صز ومزیجصهعو ۱۱0۸ بمص 0۵1 مص مت هیهت 


نتیجه گیری 

از مصادیق بارز تغييرقليم افزایش دما و کاهش بارش در اقلیم 
آینده است. بنابراین, هرگونه تغییر در وضعیت آب و هوای جهان, به 
طور مستقیم بر تولید محصولات کش‌ورزی و در نتیجه بر موضوع 
تامین امنیت غذایی در سطح جهان تاثیر گذار است. در همین راستا در 
این مطلعه به ازیابی تاثیرات تنييرقليم بر عملکرد و طول مراحل 
فنولوژیکی سیب‌زمینی در سه ایستگاه جیرفت» کهنوج و منوجان, از 
مراکز اصلی تولید سیب‌زمینی درکشور محسوب می‌شوند در آینده 
پرداخته شد. مناطق نام برده در این مطالعه جهت شبیه‌سازی اقلیم 
آینده» از داده‌های مدل گردش کلی 020۳5[۷2) تحت سناریوهای 
انتشار 10۳2.6 ,10۳4.5 و 0۳8.5 برای دوره ۲۰۱۱-۲۱۰۰ 
استفاده گردید. برای ریزگردانی خروجی مدل 020۳51۷۲2 از نرم‌افزار 
آماری 5195/1 و برای شبیه‌سازی عملکرد سیب‌زمینی از مدل رشد و 
نمو ۷۷۵۳۵۹ استفاده شد. نتایج حاصل از شبیه‌سازی پارامترهای 
دما (دمای حداقل و حداکنر» حاکی از افزايش میانگین دمای حداقل و 
دمای حداکثر به‌ترتیب از ۰/۳۸ تا ۱ درجه سانتیگراد و ۰/۳۶ ۱ ۰/۹ 
درجه سانتیگراد بود. در ایستگاه منوجان در دوره آخر روند افزایش دما 
او ستارییی 0۳45 یه ستازیوی 5 08 عکس شده است:ه 
میانگین دمای حداقل و حداکثر هر دو در این دوره کاهش می‌یابد و 
این موضوع می‌تواند نشان از بیابانی شدن هر چه بیشتر اقلیم منطقه 
در آینده دور باشد. میانگین کاهش بارش در مناطق مورد مطالعه از 


۱ میلی‌متر تا ۲۲ میلیمتر محاسبه گردید. 

نتایج کالیبراسیون و تعیین اعتبار مدل ۷۷۵۳0۵5۲ حاکی از دقت 
بالای مدل در شبیه‌سازی عملکرد و مراحل فنولوژیکی سیب‌زمینی در 
مناطق مورد مطالعه بود بطوریکه مقدار ۷18۳ برای هر دو مورد 
در مناطق تحت مطالعه زیر ۱۰ درصد بود. طول مراحل فنولوژیکی 
(کاشت تا سبزشدن» سبزشدن تا گلدهی و سبزشدن تا رسیدگی 
فیزیولوژیکی(برداشت) در دوره آینده تحت سناریوهای 136۳ کاهش 
می‌یابد. 

نتایج حاصل از شبیه‌سازی عملکرد غده سیب‌زمینی تحت شرایط 
اقلیم آینده حاکی از کاهش میزان عملکرد در مناطق مورد مطالعه بود. 
بیشترین کاهش عملکرد در دوره سوم و تحت سناریوی 6۳8.5 
مشاهده شد. بطور متوسط در مناطق مورد مطالعه به‌میزان ۱۴۷۳/۹ 
کیلوگرم در هکتار کاهش عملکرد غده سیب‌زمینی مشاهده گردید. 
کاهش طول دوره رسیدگی سیب‌زمینی بیشترین سهم را در کاهش 
عملکرد غده سیب‌زمینی دارد. نتایج بدست آمده از شبیه‌سازی عملکرد 
غده سیب‌زمینی نشان می‌دهد که به‌موازات افزايش دما میزان عملکرد 
غده نیز کاهش می‌یابد و به‌ازای هر یک درجه سانتیگراد افزایش 
میانگین دمای سالانه. عملکرد غده به‌میزان ۴/۹۹ درصد کاهش 


خواهد یافت. 


۰ نشریه بوم شناسی کشاورزی. جلد ۱۶ شماره!. بهار ۱۴۰۳ 
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